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In this contribution, we review some works related with the extraction 
of the symmetry energy parameters from isovector nuclear excitations, like 
the giant resonances. Then, we move to the general issue of how to assess 
whether correlations between a parameter of the nuclear equation of state 
and a nuclear observable are robust or not. To this aim, we introduce the 
covariance analysis and we discuss some counter-intuitive, yet enlightening, 
results from it. 


1. Introduction 

Among the widely debated questions in nuclear physics, we can mention 
the one related to the nuclear equation of state and its extrapolation to 
extreme conditions. Even restricting to zero temperature, the nuclear com¬ 
munity is still striving to determine the behaviour of the energy per particle 
in uniform matter as a function of density and neutron-proton asymmetry. 
The energy as a function of the neutron-proton asymmetry or, in turn, the 
so-called symmetry energy (cf. below for a precise defintion) is of particu¬ 
lar interest because of its impact on the physics of exotic, neutron-rich or 
neutron-deficient, nuclei. It also affects in an important fashion the proper¬ 
ties of some astrophysical compact objects like the neutron stars. Topical 
conferences address the problem of the determination of the nuclear sym¬ 
metry energy, and the reader can consult a recent topical volume to check 
the present status of our understanding |jLj. 

One can expect that the study of isovector modes of finite nuclei can shed 
light on the problem of the determination of the symmetry energy and of its 
density dependence. In the isovector collective motion, protons are displaced 
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with respect to neutrons: in other words, one creates locally a neutron- 
proton imbalance. Thus, the response of the nucleus to this perturbation is 
related to the variation of the energy density as a function of asymmetry. 
However, the nucleus is different from uniform matter: not only shell effects, 
or pairing, are expected to play a role, but isospin is not a good quantum 
number and the separation of isoscalar and isovector motion raises some 
concern. We briefly discuss some works that are related with the constraints 
on the symmetry energy emerging from the analysis of isovector properties 
of finite nuclei, and we stress the consistency of the results we have obtained. 

Then, we introduce covariance analysis as a rigorous way to determine 
whether extracting constraints on the symmetry energy from isovector ob¬ 
servables is justified or not. Until recently, covariance analysis has not been 
object of much interest by nuclear theorists. We introduce its basic concepts 
and try to make the reader familiar with the idea that there is a quanti¬ 
tative fashion to determine whether two quantities A and B are correlated 
or not, whenever calculated in a given framework. In the present context, 
A and B could be, respectively, a parameter characterizing the density be¬ 
haviour of the symmetry energy and an isovector observable. Our discussion 
will be nevertheless more general. Our analysis is based on Energy Density 
Functional (EDF) calculations [2j. We will show that the correlations that 
emerge are not universal but will, to some extent, depend on the chosen 
model. 

2. General definitons: the nuclear equation of state and the 


symmetry energy 


We assume that the nuclear systems, as any Fermi system, can be de¬ 
scribed in terms of a local energy functional. This implies that we can write 
their total energy as 



( 1 ) 


where £ is the energy density and p n and p p are, respectively, neutron and 
proton densities. Instead of p n and p p , one can use the total density p and 
the local neutron-proton asymmetry, 


0 _ Pn Pp 

P 


(2) 


In asymmetric matter, we can make a Taylor expansion of £(p, fd) in fd and 
retain only the quadratic term (odd powers of f3 are forbidden due to isospin 
symmetry), 


£{p,/3) « £o(p,P = 0) + £ sym {p)/d 2 
= £ 0 (p,fd = 0) + pS(p)p 2 . 


( 3 ) 



sym printed on April 9, 2015 


3 


The first term on the r.h.s. is the energy density of symmetric nuclear 
matter while the second term defines the the symmetry energy S(p). The 
expansion ([3]) should continue with a quartic and possibly with higher order 
terms in /?; however, the coefficient of the term in /3 4 is, to the best of our 
knowledge, found to be negligible in most models at the densities of interest 
for this work M. 

If we focus on the density dependence of S close to the usual nuclear 
density we can define 


J = S(po), 

L = 3p 0 S'(po), (4) 

where po is the saturation density for symmetric nuclear matter, po ~ 0.16 
fm~ 3 . L is often referred to as the “slope parameter”. 

We shall, in what follows, discuss how these parameters can be ex¬ 
tracted from the comparison between theoretical calculations and experi¬ 
mental measurements of the properties of isovector states. The theoretical 
calculations are done using the Hartree-Fock (HF) plus Random Phase Ap¬ 
proximation (RPA) framework. We do not provide here information about 
this well-known scheme. The details of our implementation, together with 
a general introduction, can be found in Ref. [5]. 


3. Symmetry energy extracted from giant dipole, giant 
quadrupole and pygmy dipole resonances 

3.1. IVGDR 

The case of the most collective and well known isovector giant resonance, 
namely the isovector giant dipole resonance (IVGDR), has been studied in 
Ref. [1|. In that work, the starting point is the hydrodynamical model pro¬ 
posed by E. Lipparini and S. Stringari [6j. If one denotes by rrik the k-th 
moment associated with the strength function of an external operator F, 
the IVGDR energy can be evaluated within the framework of the hydrodi- 
namical model if defined as E—\ = y/mi/m-i. The result is 


E, = 


/ 1711 

3 n 2 

^vol 

1 m -1 \ 

m(r 2 ) 

( x + 


(1 + k), 


(5) 


where 6 V oi and ^surf are the volume and surface coefficients of the macro¬ 
scopic symmetry energy and n is the well-known “enhancement factor”, 
which in the case of Skyrme forces is associated with their velocity depen¬ 
dence f2j. The coefficient 6 vo i can be identified with S(po) = J; if the nucleus 
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had a sharp surface this would be the only quantity appearing in the pre¬ 
vious expression. The nuclear surface does manifest itself in the correction 

(l + . To connect this with the microscopic symmetry energy 

5, one can rewrite this correction and assume that the r.h.s. of Eq. ([5]) in 
a heavy nucleus does not scale as y / 'S(po), but rather as \/S(p) where p is 
some value of density below the saturation density p$. In Ref. [4j it has been 
found that such a correlation between F_i (calculated within HF-RPA) and 
y/S(p) exists, with p around 0.1 fm -3 for heavy nuclei (in agreement with 
[7]). By exploiting this correlation and inserting the experimental value for 
the IVGDR energy in 208 Pb, it has been found that 

23.3 MeV < 5(0.1) < 24.9 MeV. (6) 


3.2. IVGQR 

It is a natural question to ask whether the isovector giant quadrupole 
resonance (IVGQR) provides a consistent extraction of the density depen¬ 
dence of the symmetry energy. Until recently, the experimental properties 
of the IVGQR had not been determined very accurately; this goal has been 
achieved using a very intense and polarized photon beam at the HlyS facil¬ 
ity IS] - In Ref. [9], a comprehensive theoretical analysis has been performed, 
based both on RPA calculations and on a macroscopic interpretation. The 
energy of the IVGQR receives contribution from unperturbed particle-hole 
(p-h) configurations at 2hu excitation energy, plus some correlation energy 
related to the isovector residual interaction. This idea has been implemented 
in Ref. [9], with mild assumptions and taking care of the fact that the un¬ 
perturbed p-h energy can be related to the effective mass and, in turn, to 
the isoscalar GQR energy. The main result is 


-^TVGQR ~ 2 


(£isgqr ) 2 0 (3S(p) 

2 + A 2 / 3 U Foo 



(7) 


where is the Fermi energy for symmetric nuclear matter at saturation 
density, and S(p) is the symmetry energy at some average nuclear density. 
If we take for this the same value that we have adopted in the above discus¬ 
sion for the IVGDR, that is, 0.1 fm -3 , we can reproduce the experimental 
IVGQR energy; to turn it around, from the two experimental IVGDR and 
IVGQR energies we can derive consistent values for the value of 5(0.1). 

In Ref. [9] it has been checked, in addition, that both Skyrrne and 
relativistic mean-field models follow quite well the scaling predicted by Eq. 
Q. In fact, new Skyrrne interactions have been fitted in Ref. [9] using the 
same protocol as in Ref. |10], where the new set SAMi has been introduced; 
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moreover, new effective Lagrangians have been fitted along the line of the 
DD-ME one HU. These sets are characterized by different values of J and of 
the effective mass m*/m. The microscopic results obtained with sets having 
the same value of mn*/m have been used as follows. The IVGQR energy can 
be reproduced only with forces having a specific combination of J and T; if 
we assume a value of J = 32±1 MeV we extract 

L = 37 ± 18 MeV; (8) 

Ai?( 208 Pb) = 0.14 ±0.03 fin. (9) 

The second line corresponds to the neutron skin in 208 Pb, that is well known 
to be correlated with L jT2j . 


3.3. PDR 

Among isovector modes, the so-called “Pygmy Dipole Resonance” (PDR) 
has captured noticeable interest in the last decade. This definition is not 
free from ambiguities if used to label generically the dipole strength below 
the IVGDR. Experimentally, low-lying dipole strength has been found in 
several nuclei 113111-41 . Typically, the PDR strength may arrive up to a few 
% of the dipole EWSR. 

In Ref. [T5] , it has been proposed that the fraction of EWSR exhausted 
by the PDR and the slope parameter L defined in Eq. <s> are correlated. 
This correlation could be explained if PDR is a mode related to the oscilla¬ 
tion of the excess neutrons, whose dynamics is decoupled from the IVGDR 
(see Ref. [16] for a transparent interpretation). However, this picture may 
break in some cases. The collectivity of the PDR seems to be somewhat 
model-dependent, and the states in that energy region have also a mixed 
isovector/isoscalar character [U, jT8] 19]. A toroidal component has been 
found in the calculations of Refs. 120 EH- Despite these well-taken warn¬ 
ings, the values of the slope parameter and neutron skin (for 208 Pb) that 
have been extracted in UU, namely 

L = 64.8 ± 15.7 MeV; 

A R ( 208 Pb) = 0.194 ± 0.024 fm, 

are consistent with the previous values given in Eq. Q and Q 

3-4- Total dipole polarizability 

The total dipole polarizabilty ap> is proportional to the inverse-energy 
weighted sum rule m_i of the dipole operator, the exact relationship being 

87re 2 


( 10 ) 

( 11 ) 


OtD = 


9 


rri-\. 
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The correlations between this quantity, calculated by means of a large bunch 
of EDFs, and either the slope parameter L or the neutron skin has been 
analyzed in Ref. [22] • In that work, however, the correlation shows up very 
clearly only within families made up with similar models. In fact, with the 
help of the droplet model (DM), it has been shown in Ref. [23] that the 
slope parameter L is well correlated with the product of J times the dipole 
polarizability. The formula suggested by the DM is 


DM 

a D 


vre 2 A(r 2 ) 
J 


1 + 


_ |\„DM I /T e 2 Z A „surf 
5^'np T- Y 5 70J np 

2 (r 2 )i/ 2 (/ — I c ) 


(13) 


where the quantities appearing in the second term within square brackets 
are defined in [;23j and shown to not vary appreciably among models (at 
least in heavy systems like 208 Pb [24]). Microscopic calculations obey such 
kind of scaling quite well and allow extracting, assuming J = 31 ± 2 MeV, 
the value 

L = 43 ± 16 MeV. (14) 


4. Covariance analysis 

All the above discussion is simply based on the empirical appearance 
of a linear correlation between two observables when they are calculated 
using several models like EDF parameterizations. This kind of analysis is 
not based on any statistical assumption (models are assumed to be inde¬ 
pendent) but could be justified if the correlation under study is suggested 
by macroscopic models. A different strategy to judge about correlations be¬ 
tween observables is to use covariance analysis. The advantage is that there 
is a more rigorous statistical foundation for such a method; however, this 
can be used only to judge whether correlations exist within the framework of 
a given model, whose parameters are varied without changing the ansatz or 
the fitting protocol of such parameters. Although this will not be our focus 
here, covariance analysis is interesting for several more reasons. It allows 
estimating the theoretical errors on the model parameters and thus, decide 
if some of them is underconstrained or even redundant. 

For the sake of brevity we do not give here a thorough explanation of 
the method of covariance analysis. We refer to our recent work [25] for 
better explanations, reminding also that many textbooks contain a more 
exhaustive treatment of the formalism (see e.g. [26]). and that an excellent 
introduction for nuclear theorists is given in Ref. m- 
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4-1■ Summary of relevant formulas 

Let us consider a model, like an EDF one, characterized by n parameters 
p = (pi, ...,p n ). Observables (.A) are functions of these parameters. Usually 
one builds the optimal model, in the space defined by the n parameters, 
by x 2 minimization. Let us also assume that the x 2 is a we ll behaved, 
analytical hyper-function of the parameters around their optimal value po, 
and that it can be approximated there by a Taylor expansion, namely 

1 n 

X 2 (P) - X 2 (Po) ~ x ^{Pi ~ P 0 t)d Pl d Pj x 2 (Pj ~ Poj) 

hJ 

n 

= ~ P^) (Pj ~ Poj) , (15) 

hj 


where we have defined the curvature matrix £~ 1 and its inverse which is the 
covariance (or error) matrix £. 

Let us now expand an observable A(p) around the minimum po assuming 
a smooth behavior and, therefore, neglecting second order derivatives: 

A(p) = ^(p) + (P - Po)<9p-4(p) |p=p 0 = A 0 + (p — Po)A 0 . (16) 

The covariance between two observables is defined as 


Cab = (A(p) - A)(B(p) - B) ■ (17) 


Using the above expansion for the observables, and the fact that the param¬ 
eters they depend upon are expected to have a Gaussian distribution with 
curvature matrix T _1 , one can write 


Cab ~ 


&4(P) 

dpi 






P=PO 


dB( p) 


dpi 


P=P0 


(18) 


The variance of A is, then, given by Caa ■ One may also calculate the 
Pearson-product moment correlation coefficient between those observables, 
that is, 

_ Cab 
cab 

V^AA^BB 

cab = 1 means complete correlation between observables A and B, whereas 
— 1 means complete anti-correlation and cab = 0 means no correlation at 
all. 
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(a) 



(b) 



Fig. 1. Values of the Pearson-product correlation coefficient calculated among pairs 
of observables. Panels (a) and (b) corresponds to two different kinds of mimimiza- 
tion procedures, that are discussed in the text. 


4-2. Some results from the covariance analysis 

Several results obtained via the covariance analysis have been presented 
in Ref. [25]. Here, we only focus on some illustrative findings that are of 
relevance in connection with our previous discussion on symmetry energy 
and other isovector observable. 

We have fitted a functional called SLy5-min, that has been built so that 
to be as close as possible to the original SLy5 force [28] [29] . Details and 
differences are highlighted in Ref. [25]. We now illustrate what happens to 
correlations between observables when one varies the fitting protocol. 

The original SLy5 functional, and the present Sly5-min functional as 
well, has been fitted by including in the x 2 as pseudo-data a set of values 
for the neutron matter equation of state derived by some ab-initio calcula¬ 
tion with realistic forces [30] . In the variant that we label as SLy5-a, we have 
changed in the x 2 the weight associated with these pseudo-data, that is, with 
the equation of state of neutron matter. In particular, We have increased 
the value of the error on these points, A^(p, 5 = 1), from 0.1 x ^(p, <5 = 1) - 
that corresponds to a 10% relative error - to 0.5 x ® (p, 6 = 1). The Pearson- 
product correlation coefficients of this fit are shown in panel (a) of Fig. [l] 
One can notice a strong correlation between the isovector observables, that 
are, the neutron radius of 208 Pb, S'(po), L and m_i(IVGQR). This result 
clearly indicates that when a constraint on a property A is relaxed, corre¬ 
lations with other related observables B not included in the fitting protocol 
become larger. 

The second variant we have built is denoted as SLy5-b. In this case, we 
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have kept all terms in the y 2 as i n SLy5-min except the equation of state 
of neutron matter that now is not included at all, and we added instead a 
very tight constraint on the neutron skin thickness of 208 Pb: we have chosen 
Ar np = 0.160T0.001 fm. Then, panel (b) of Fig. [I] shows another interesting 
outcome: A r np displays an almost vanishing correlation with all the other 
quantities. This indicates that when a property A is tightly constrained in 
the fitting protocol, correlations with other observables B become very small. 


5. Conclusions 

The most important outcomes of our recent works on the constraints on 
the density dependence of the symmetry energy imposed by the properties 
of isovector excitations are 


if one looks at correlations between observables and values or deriva¬ 
tives of the symmetry energy, calculated using a large set of different 
EDFs, consistent values of L and of the neutron skin thickness in 
208 Pb are extracted from the study of the IVGDR, IVGQR, PDR and 
total dipole polarizability [cf., in particular, Eqs. ([8]), (10), (14) as 
well as Eqs. <©>©]; 


• some of these correlations are suggested by macroscopic models and, in 
general, we believe that it would be desirable that correlations emerge 
from a sound physical picture rather than simply from a numerical 
analysis; 


• models may, or may not, display the correlations that one expects on 
physical grounds because of the strong impact of the fitting protocols. 
We have shown that correlations of an observable A with related ob¬ 
servables B may emerge (vanish) if the constraint put on A in the 
fitting protocol is negligible (very tight). 
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